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Abstract 



The prospects at the CERN LHC are investigated for the discovery via 
decays into muon pairs of neutral Higgs bosons in the minimal supergravity 
model. Promising results are found for the CP-odd pseudoscalar (A^) and 
the heavier CP-even scalar (H^) Higgs bosons at large tan/? = v^jvx. For 
tan/3 ^ 10 and 100 GeV ^ mA,mH ^ 550 GeV, the A^ and the masses 
could be precisely reconstructed from the dimuon invariant mass. 
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Typeset using REVT^ 



I. INTRODUCTION 



In the Standard Model (SM) of electroweak interactions, the masses of gauge bosons and 
fermions are generated by a scalar field doublet. After spontaneous symmetry breaking, a 
neutral CP-even Higgs boson (hg^) remains as a physical particle. The scalar sector of the 
SM is unstable to radiative corrections that arise when the SM is coupled to ultraheavy 
degrees of freedom such as those in a grand unified theory (GUT). Therefore, there is great 
interest in extensions of the SM that can solve this problem. 

A supersymmetry (SUSY) between fermions and bosons provides a natural explanation 
of the Higgs mechanism for electroweak symmetry breaking (EWSB) in the framework of 
a grand unified theory (GUT) and preserves the elementary nature of Higgs bosons. For 
the particle content of the minimal supersymmetric standard model (MSSM) [|I|, the evolu- 
tion of gauge couplings by renormalization group equations (RGBs) is consistent with a 
grand unified scale at Mqut ~ 2 x 10^^ GeV and an effective SUSY mass scale in the range 
Mz < MsusY ~ 10 TeV 0. With a large top quark Yukawa coupling (Yt) to a Higgs boson 
at the GUT scale, radiative corrections drive the corresponding Higgs boson mass squared 
parameter negative, spontaneously breaking the electroweak symmetry and naturally ex- 
plaining the origin of the electroweak scale. In the minimal supersymmetric GUT with a 
large Yt, there is an infrared fixed point (IRFP) at low tan/3; the top quark mass is 
correspondingly predicted to be = (200 GeV) sin/? [Q, and thus tan/5 ~ 1.8 for nit = 175 
GeV. At high tan/?, another quasi-IRFP solution {dVt/dt ~ 0) exists at tan/? ~ 56. 

The Higgs sector of a supersymmetric theory must contain at least two SU (2) doublets 
for anomaly cancellation. In the minimal supersymmetric standard model (MSSM), the 
Higgs sector has two doublets (pi and 02 that couple to the ts = —1/2 and = +1/2 
fermions, respectively. After spontaneous symmetry breaking, there remain five physical 
Higgs bosons: a pair of singly charged Higgs bosons H^, two neutral CP-even scalars 
(heavier) and (lighter), and a neutral CP-odd pseudoscalar A'^. The Higgs potential is 
constrained by supersymmetry such that all tree-level Higgs boson masses and couplings 
are determined by just two independent parameters, commonly chosen to be the mass of 
the CP-odd pseudoscalar (m^) and the ratio of vacuum expectation values (VEVs) of Higgs 
fields (tan/? = V2/V1). 

Extensive studies have been made of the prospects for the detection of MSSM Higgs 
bosons at the LHC ||7|-[T3|. Most studies have focused on the SM decay modes — ^ 77 



H^, or v4°) and (p ZZ or ZZ* 4/ (</> = or /i°) that are detectable above 
background. For tan (3 close to one, the detection modes —>■ Zh^ —>■ l^l^bb or l^l~Tf |1^ 
and —>■ h^h^ '-f-fbb [O may provide promising channels to simultaneously discover 



two Higgs bosons of the MSSM. There are regions of parameter space where rates for Higgs 
boson decays to SUSY particles are large and dominant. While these decays reduce the 
rates for the standard modes, making conventional detection of Higgs bosons more difficult, 
they also open up a number of new promising modes for Higgs detection |10| . 



For large tan/3, the rf decay mode [^|TT],|T^ is a promising discovery channel for the 
A^ and the in the MSSM. It was also suggested that neutral Higgs bosons might be 
observable via their bb decays [Il5| , |l6| in a large region of the (mA,tan/3) plane, provided 
that sufficient 6-tagging capability can be achieved. However, simulations with ATLAS 



detector performance concluded that detection of the bb channel would be difficult |17 
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Recently, the muon pair decay mode was proposed [p^pl]| and confirmed [T^ to be a 
promising discovery channel for the and the H^. Although the fifi channel has a small 
branching fraction, this is compensated by the much better achievable mass resolution with 
muon pairs. For large tan/3, the muon pair discovery mode might be the only channel at 
the LHC that allows precise reconstruction of the and the masses. 

In supergravity (SUGRA) models [0, supersymmetry is broken in a hidden sector with 
SUSY breaking communicated to the observable sector through gravitational interactions, 
leading naturally but not necessarily pO| to a common scalar mass (mo), a common gaugino 
mass (mi/2), a common trilinear coupling (Aq) and a common bilinear coupling (Bq) at the 
GUT scale. Through minimization of the Higgs potential, the B coupling parameter of the 
superpotential and the magnitude of the Higgs mixing parameter /i are related to the ratio 
of VEVs of Higgs fields (tan/3 = V2/V1) and to the mass of the Z boson (M^). The SUSY 
particle masses and couplings at the weak scale can be predicted by the evolution of RGBs 
from the unification scale P,PT|. 

For tan/3 close to 1.8, the lower IRFP, masses of the CP-odd pseudoscalar (m^) and 
the heavier CP-even scalar (m^/) are large over most of the minimal supergravity parameter 
space. Then the couplings of the lighter scalar are similar to those of the SM Higgs boson; 
the could be the only neutral Higgs boson observable at the CERN LHC. For large tan (3, 
the A^ and the can become light |2^J2H] and potentially visible at the LHC. 



Recent measurements of the b ^ S'j decay rate by the CLEO p4[ and LEP collaborations 
25| place constraints on the parameter space of the minimal supergravity model p6[. It 



was found that 5 — 57 excludes most of the minimal supergravity (mSUGRA) parameter 
space when tan/3 is large and /i > [^]. Therefore, we will choose /i < in our analysis. 
However, our results are almost independent of the sign of /z. 

The remainder of this article addresses the prospects of discovering the neutral Higgs 
bosons in the mSUGRA model via their decays into muon pairs at the LHC. The cross section 
and branching fractions are presented in Section II. The observability of the dimuon signals 
is discussed in Section HI for conservative assumption about the detector mass resolution. 
Promising conclusions are drawn in Section IV. 



II. CROSS SECTION AND BRANCHING FRACTION 

We evaluate SUSY mass spectra and couplings in the minimal supergravity model with 
four parameters: mo, mi/2, Aq and tan /?, and with the sign of the Higgs mixing parameter /x. 
Since Aq mainly affects the masses of third generation sfermions, which do not significantly 
affect our analysis, we take = in our calculations. 

The mass matrix of the charginos in the weak eigenstates {W^, H"^) has the following 
form m 

^( M2 V2Mwsm(3\ 
^ \V2Mw cosp -n J ' ^ ^ 

This mass matrix is not symmetric and must be diagonalized by two matrices [||] . The sign 
of the /i contribution in Eq. (|^) establishes our sign convention for fi. 

We calculate the masses and couplings in the Higgs sector with one-loop corrections 
from both the top and the bottom Yukawa interactions in the RGE-improved one-loop 
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effective potential at the scale Q = ^Jfnq^^rni^. With this scale choice, the numerical 
value of the CP-odd Higgs boson mass (m^) at large tan/3 |p2| , p3| is relatively insensitive 
to the exact scale choice and the loop corrections to are small compared to the tree 
level contribution. In addition, when this high scale is used, the RGE improved one-loop 
corrections approximately reproduce the dominant two-loop perturbative calculation of the 
mass of the lighter CP-even Higgs scalar (m/j). Our numerical values of nih are very close to 
the results of Ref. where somewhat different scales higher than Mz have been adopted 
in evaluating the effective potential. 

The cross section of pp —>■ (p ^ fip, + X {(p = A^, if °, or h^) is calculated for the two 
dominant subprocesses gg (p and gg —>■ (pbb, and multiplied by the branching fraction 
of the Higgs decay into muon pairs B{(f) ^p,). The parton distribution functions of 
CTEQ3L ||2^ are used to evaluate the pp ^ (p + X cross section with A4 = 0.177 GeV and 
Q = Mgg = the invariant mass of the gluons. We take Mz = 91.19 GeV, sin^ 9w = 0.2319, 
M]y = Mz COS 9w, mh{pole) = 4.8 GeV, and rritipole) = 175 GeV. 

At the LHC energy, the SM Higgs boson is produced dominantly from gluon fusion; 
vector boson fusion is also relevant if the Higgs boson is heavy. In the MSSM, gluon fusion 
{gg — *• (p) is the major source of neutral Higgs bosons for tan j3 less than about 4. If tan /3 
is larger than about 10, neutral Higgs bosons in the MSSM are dominantly produced from 
6-quark fusion {bb cp) [3C]. The cross section of gg — > (phb is a good approximation to 
the 'exact' cross section |^ of bb —>■ (p for M^ less than about 500 GeV. Since the Yukawa 
couplings of (pbb are enhanced by 1/ cos/5, the production rate of neutral Higgs bosons is 
usually enhanced at large tan /?. For rriA larger than about 150 GeV, the couplings of the 
lighter scalar to gauge bosons and fermions become similar to those of the SM Higgs 
boson. Then gluon fusion is the major source of the even if tan/3 is large. 

The QCD radiative corrections to the subprocess gg cp are substantial pT|,|52|; the 



corrections to gg (pbb are still to be evaluated. To be conservative, we take a K-factor 
of 1.5 for the contribution from gg ^ (p and a K-factor of 1.0 for the contribution from 



gg — > (pbb. For the dominant Drell-Yan background |18,11,13|, we adopt the well known 



K-factor from reference 33 



With QCD radiative corrections to (p ^ bb ||3J], the branching fraction of ^ /i/i is 
about m^/SmKmb) ~ 2 x 10~^ when the bb mode dominates Higgs decays, where 3 is a 
color factor of the quarks and mb{mb) is the running mass at the scale mb. The branching 
fraction of — >■ /i/2 is always about 2 x 10~^. The branching fractions for fifi and 

— > up, are always in the range 1.5 — 2.5 x 10"^ when tan/3 ^ 10, even when A^ and 
can decay into SUSY particles. For ttia less than about 80 GeV, the decays dominantly 
into h'^h^, AO AO and ZA°. 

In Figures 1(a) and 1(b), we present the cross section of the MSSM Higgs bosons at the 
LHC, pp (p fip + X, as a function of ttia for rriq = irtg = — = 1 TeV, (a) tan/? = 15 
and (b) tan/3 = 40. As tan/3 increases, the production cross section is enhanced because for 
tan P > 10, the production cross section is dominated by gg (pbb and enhanced by the (pbb 
Yukawa coupling. Also shown is the same cross section for the SM Higgs boson /i^j^j with 
mhgj^j = rriA- For nihgj^j > 140 GeV, the SM hgj^^ mainly decays into gauge bosons, and the 
branching fraction B{h^gj^,^ fip) drops sharply in this mass region. 

We present the cross section of pp ^ (p ^ fifi + X as a function of tan /? for /i < 0, 
mo = 300 GeV and various values of mi/2 in Figure 2. The regions in Fig. 2 with dark 
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shading denote the parts of the parameter space that do not have the hghtest neutrahno 
(Xi) as the hghtest supersymmetric particle (LSP). As tan/? increases, the production cross 
section is enhanced because (i) is reduced; and (ii) for tan /? > 10, the production cross 
section is enhanced by the (pbb Yukawa couphng. Also shown in Fig. 2 are the pseudoscalar 
mass (m^) versus tan (3 and the invariant mass distribution {da/dM^fi) of the background 
from the Drell-Yan process at M^^ = itla- 

In Figure 3, we show the cross section of — > </> — > /x/i + X in fb, as a function of mi/2, 
with = 14 TeV, for > 0, mo = 300 GeV, and (a) tan/5 = 1.8, (b) tan^ = 10, (c) 
tan/? = 35, and (d) tan/? = 50. Also shown are irtA versus mi/2 and the invariant mass 
distribution {da /dM^^p) of the background from the Drell-Yan process at M^^ = rriA. The 
regions in Fig. 3 with dark shading denote the parts of the parameter space that do not 
satisfy the following theoretical requirements: tachyon free, and the lightest neutralino as 
the LSP. For tan/? ~ 1.8, rriA is very large and the cross section of the A^,H^ signal is 
much smaller than the background. For tan/5 ^ 35, the cross section of the A^,H^ signal 
is greatly enhanced and can become slightly larger than the background. For tan/5 ~ 50 
and mi/2 < 200 GeV, can become close to mz and a peak appears in the invariant mass 
distribution of the background from the Drell-Yan process. 



III. OBSERVABILITY AT THE LHC 



We define the signal to be observable if the Na lower limit on the signal plus background 
is larger than the corresponding upper limit on the background [|7|,^, namely. 



L{a, + ab) - NJL{a 



> L(Tb + N\ La, 



(2) 



which corresponds to 



a. > 



Ar2 

T 



1 + 2JLab/N 



(3) 



Here L is the integrated luminosity, as is the cross section of Higgs signal, and ab is the back- 
ground cross section within a bin of width ±AM^^ centered at Mf, N = 2.32 corresponds 
to a 99% confidence level and = 2.5 corresponds to a 5a signal. We take the integrated 
luminosity L to be 300 fb^^ 

To study the observability of the muon discovery mode, we consider the background from 
the Drell-Yan (DY) process, qq —>■ Z,'~f ^ njl, which is the dominant background. We take 
AM^fi to be the larger of the ATLAS muon mass resolution (about 2% of the Higgs bosons 
mass) |]12|,|13| or the Higgs boson width. The CMS mass resolution will be better than 2% 
of for ^ 500 GeV [|ll|,0. Therefore, the observability will be better for the CMS 
detector. The minimal cuts applied are (1) Pt(/^) > 20 GeV and (2) \f]{fi)\ < 2.5 for both 
the signal and background. For ^ 130 GeV, uia and ttih are almost degenerate while 
for niA ^ 100 GeV niA and nih are very close to each other |T8|JTT[] . Therefore, we add up 
the cross sections of the and the for ttia < 100 GeV and those of the A^ and the 
for ruA > 100 GeV. 

The 5a discovery contours for the MSSM Higgs bosons at y/s = 14 TeV with an integrated 



luminosity L = 300 fb are shown in Figs. 1(c) and 1(d) for rriq 



-fi = 1 TeV 
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and rriq = rrig = —jj, = 300 GeV. The discovery region of /i/i is slightly enlarged for 

a smaller uia, but the observable region of nJi is slightly reduced because the lighter 
top squarks make the and the h*^ lighter; also the H%h coupling is enhanced while the 
h^bb coupling is reduced. 

In Figure 4, we present the LHC discovery contours in the minimal supergravity model, 
for (a) the mi/2 versus tan/3 plane with mo = 150 GeV, (b) the mi/2 versus tan/3 plane with 
mo = 500 GeV, (c) the mi/2 versus mo plane with tan/? = 15, and (d) the mi/2 versus mo 
plane with tan /5 = 40. The discovery region is the part of the parameter space between the 
curve of square symbol and the dash line. The QCD radiative corrections to background 
from the Drell-Yan process are included. The regions in Fig. 4 with dark shading denote 
the parts of the parameter space that do not satisfy the following theoretical requirements: 
electroweak symmetry breaking (EWSB), tachyon free, and the lightest neutralino as the 
LSP. Also shown are the mass contours for uia = 100 GeV (dashed), 500 GeV (sohd) and 
1000 GeV (dot-dashed). For tan/? > 40, it should be possible to observe the AP or the 
with a mass up to about 550 GeV. 

There are a couple of interesting points to note: (i) an increase in tan (3 leads to a larger 
nth but a reduction in itla and mj^; (ii) an increase in mo or in mi/2 raises m^, mu and 
masses of the other scalars significantly. 



IV. CONCLUSIONS 

The muon pair decay mode is a very promising channel to discover the neutral Higgs 
bosons of minimal supersymmetry and minimal supergravity. The discovery region of the 

mode is slightly smaller than the rf channel but /i/i allows precise reconstruction of the 
Higgs boson masses. The and should be observable in a large region of parameter 
space with tan/? > 10. The should be observable in a region with niA < 120 GeV and 
tan /? ^ 5. 

The observable regions of the parameter space are found to be 

mo = 150 GeV : mi/2 < 400 GeV and tan /? > 12 

mo = 500 GeV: mi/2 < 1 TeV and tan/? > 28 (4) 

For two specific choices of large tan /3, the observable regions are 

tan /? = 15 : mi/2 < 200 GeV and mo < 200 GeV 

tan /? = 40 : mi/2 < 600 GeV and mo < 800 GeV (5) 

All observable regions are nearly independent of the sign of /x. 

For tua ~ 200 GeV and tan/? > 25, L = 10 fb"^ would be sufficient [18| to obtain 
Higgs boson signals with a statistical significance larger than 7. For M^^ close to the Mz, 
the signal is marginal because it appears on the shoulder of the huge Z peak. Adequate 
subtraction procedures are required to extract the signal in this region. 
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FIG. 1. The total cross section of — > ^ — > fip, + X in fb at ^/s = 14 TeV, as a 

function of in the MSSM, for rrig = rriq = —fi = 1 TeV, (a) tan/? = 15 and (b) tan/3 = 40. 
Also shown are the cross section for the SM Higgs boson with rrih^j^j = tua (dotted). The 5(T 
contours at the LHC with L = 300 fb~^ are shown for (c) rrig = rriq = —fi = 1 TeV, and (d) 
rrig = rriq = —fi = 300 GeV. The discovery region is the part of the parameter space with light 
shading. 
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FIG. 2. The total cross section of ^ fifi + X in fb at ^/s = 14 TeV, as a function 

of tan/3, for /i < 0, mo = 300 GeV, and (a) my2 = 150 GeV, (b) mi/2 = 300 GeV, (c) mi/2 = 500 
GeV, and, (d) mi/2 = 800 GeV, Also shown is the invariant mass distribution [da /dMfj^p) of the 
Drell-Yan background (dashed) in fb/GeV for Af^^ = m^. The region with dark shading denotes 
the part of the parameter space that does not have the lightest neutralino (x?) as the LSP. 
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FIG. 3. The total cross section of pp A^,H^ /u/i + X in fb at ^/s = 14 TeV, as a 
function of mi/2, for /i < 0, = 300 GeV, (a) tan/? = 1.8, (b) tan/3 = 10, (c) tan/3 = 35, 
and (d) tan/3 = 50. Also shown is da/dM^fi of the Drell-Yan background (dashed) in (fb/GeV) 
for M^fi = niA- The regions with dark shading are the parts of the parameter space excluded 

by theoretical requirements. The region excluded by the m + > 85 GeV limit from the chargino 

^1 

search ^ at LEP 2 is indicated. 
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FIG. 4. The 5a contours at the LHC for an integrated luminosity (L) of 300 fb~^ in (a) the 
"^1/2 versus tan/3 plane with niQ = 150 GeV, (b) the mi/2 versus tan/3 plane with mo = 500 
GeV, (c) the mi/2 versus mo plane with tan/3 = 15, and (d) the mi/2 versus mo plane with 
tan/3 = 40. Also shown are the mass contours for niA = 100 GeV (dashed), 500 GeV (solid) 
and 1000 GeV (dot-dashed). The discovery region is the part of the parameter space with light 
shading. The regions with dark shading are the parts of the parameter space excluded by theoretical 
requirements. The region excluded by the m + > 85 GeV limit from the chargino search |35] at 



LEP 2 is indicated. 
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